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Abstract 

Heat fluctuations are studied in a dissipative system with both me- 
chanical and stochastic components for a simple model: a Brownian 
particle dragged through water by a moving potential. An extended 
stationary state fluctuation theorem is derived. For infinite time, this 
reduces to the conventional fluctuation theorem only for small fluc- 
tuations; for large fluctuations, it gives a much larger ratio of the 
probabilities of the particle to absorb rather than supply heat. This 
persists for finite times and should be observable in experiments sim- 
ilar to a recent one of Wang et al. 

PACS: 05.40.-a, 05.70.-a, 44.05. +e, 02.50.-r 

1 Introduction 

There is a lack of unifying principles in nonequilibrium statistical mechanics, 
compared to the equilibrium case. So it is not surprising that the fluctuation 
theorem has received a lot of attention, as it gives a property of fluctuations 
of entropy production for a large class of systems, possibly arbitrarily far 
from equilibrium^ |2IIHIll3inilIllEllHllIII]- 

In terms of heat rather than entropy production, the conventional sta- 
tionary state fluctuation theorem (SSFT) 1 states that the probability P T (Q T ) 

1 We call equation the "conventional SSFT" to distinguish it from the extended 
SSFT presented in this paper. 



to find a value of Q T for the amount of heat dissipated in a time interval r, 
satisfies, in a nonequilibrium stationary state [HE], 
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Pr(-Qr) 
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where ~ indicates the behavior for large r. Here (3 = (A^T) -1 , with ks 
Boltzmann's constant and T the (effective) temperature of the system. In 
contrast, the transient fluctuation theorem (TFT) considers fluctuations Q T 
in time, when the system is initially in equilibrium 2J. These theorems were 
first demonstrated in deterministic many particle systems in an external field 
[HElEj but later in stochastic systems as welljUE]. 

In this context, it is customary to identify Q T /T as a (generalized) en- 
tropy production. Then, equation (JTJ) is interpreted as a theorem for entropy 
production arbitrarily far from equilibrium. Furthermore, the fluctuation 
theorem holds for arbitrary values of Q T , i.e., also far from its average. Hence, 
it is often referred to large deviation theorem. 

Recently, a laboratory experiment was carried out by Wang et al. jHj. They 
measured fluctuations in the work done on a system in a transient state of 
a Brownian particle in water, subject to a moving, confining potential. The 
TFT was confirmed. 

In the deterministic models, dissipation is often modeled by including 
a damping term in the equations of motion, chosen such as to keep, e.g., 
the total energy of the system fixed. This damping term is a mechanical 
expression for what in reality is an external thermostat. Any work done by 
the external field on the system is absorbed by this "internal" thermostat. 
So the external work done on the system equals the heat dissipated in the 
system. 

While Wang et al. intended to study the entropy production (or heat) 
fluctuations, in fact, the work fluctuations were studiedJTU]. In contrast to 
the above sketched purely deterministic models, work fluctuations differ from 
heat fluctuations in their system due to the presence of a confining potential. 
Thus, some of the external work done is converted into potential energy 
and only the rest is converted into heat. In fact, in this system the work 
fluctuations in the stationary state satisfy the conventional SSFTjTQJ ITT] , 
but the heat fluctuations do not, as we shall show. As it turns out, having 
a deterministic together with a stochastic motion results in a behavior of 
the heat fluctuations that coincides with the conventional SSFT only for a 
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restricted set of small fluctuations and in a very different behavior for larger 
ones. 



2 Work fluctuations — conventional SSFT 

We first discuss the work-related fluctuation theorem in the experiment of 
Wang et ai, which was treated theoretically in Refs. [TUJ HH H2] ■ The Brow- 
nian particle in a fluid, subject to a harmonic potential moving with constant 
velocity v*, was described by an overdamped Langevin equation: 

£ = -(*-*:> +6. (2) 

Here, x t is the position of the particle at time t, x£ = v*t is the position 
of the minimum of the harmonic potential at time t [cf. equation (J7 7 ))]. and 
(t is a fluctuating force with zero mean and a delta function correlation in 
time. We remark that the relaxation time of the position of the particle 
has been set equal to one. Also, we set ksT = 1, so (CtCs) = 2d>(t — s). 2 
In reference JT0], it was shown that equation (j2J) is solvable in a co-moving 
frame, in which it reduces to a standard Ornstein-Uhlenbeck process. Thus, 
the stationary probability distribution and Green's function are known, and 
are both Gaussian in x t . The work is the total amount of energy put into 
the system in a time r. This is a fluctuating quantity, given by 3 



f [-(**- *t*)R (3) 
Jo 



Here, the time t = denotes the initial time of an interval of length r in the 
stationary state. W T is a linear function of the positions x t , and since those 
have a Gaussian probability distribution function, so does W T . When the 
mean and variance of the probability distribution function are computed 
[using the stationary solution and Green's function of equation (J2J)], one 

finds pn 

lim — In 



P?(PWT) ^ (4) 



2 Compared to Ref. /c^T = 1, the harmonic force constant k = 1 and the friction 
coefficient a — 1. So the energy unit is fcgT, the time unit is the relaxation time a/k, and 
the length unit is the thermal width y^ksT/k. 

3 A force on the particle implies an opposite force on U, which must be overcome by an 
external force to keep U moving. The work done by this force is W T in equation P|) 
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Here, p is a scaled value of W T , defined as p = W T / (W T ), such that (p) = 1. 
We also wrote 

(W T ) = wr, (5) 

with w the average work rate, which is independent of r in the stationary 
state. In the current units, w = |v*| 2 . Equation (£Q) is, for the work fluc- 
tuations, a more careful formulation of the SSFT in equation (JIJ). A work 
related TFT also holds [IQJ HD E21- 



3 Heat fluctuations 



3.1 Fourier transform of the distribution 

We now turn to the heat SSFT. The heat Q T is that part of the work W T 
that goes into the fluid. Some work is also stored in the potential, so 

Q T = W T - AU T , (6) 

where AU T is the change in potential energy of the particle in a time r, 

AU T = U T - U Q , (7) 

with 

U t ^ t -^\\ (7') 

This form of U% makes Q T nonlinear in x t . As a result, the probability 
distribution function P T (Q T ) of Q T need not be Gaussian. Nonetheless, it is 
possible to compute its Fourier transform. 

The Fourier transform of P T (Q T ), defined as 

/oo 
dQ T e l ^P T (Q T ), (8) 
-oo 

is computed by writing P T as [using equations © and (JTj)] 

PAQr) = JJd* rfx T Pf- X0 ^(Q r + AU T , x , x r ), (9) 

where p^' x 0' x ^ j s the joint distribution of the work W T , the positions xo 
and x T at the beginning and at the end of the time interval r, respectively. 
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This distribution is Gaussian because W T , x , and x r are all linear in x t . 
When equation is inserted into equation (|SJ), a seven dimensional Gaussian 
integral is left, which after some algebra yields 



Pr(Q) 



cxp 



{wq(i -q) r - T 
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(10) 



[l+(l_ e -2r)g2]3/2 



Once P T {q) has been transformed back, one considers 



/t(p) = — ln 



1 f P T {jpWT) 



(11) 



wt \_P T (—pwr) 



Here, p is a scaled value of Q r , defined as 



V = Qt/(Qt), 



i.e., (p) = 1. We also used that 



<Qr) = (Wr) - (AC/ r ) = »T 



by equation (JHJ), since (AU T ) = in the stationary state. The scaled loga- 
rithmic ratio / T (p) should be equal to p for r — > oo when the conventional 
SSFT holds. 

As far as we know, there is no exact result for the inverse Fourier trans- 
form of P T {q) in equation (fTU|) in terms of known functions. Therefore, a 
completely analytic treatment did not seem feasible. Instead, we used first a 
numerical method, the fast Fourier transform algorithm^!!] > to invert equa- 
tion (JlOj) . The resulting probability distribution function P T as well as the 
corresponding f T have been plotted in figure ^ These results do not agree 
very well with the straight line with slope one, which should be approached 
for large r if the conventional SSFT were to hold. One might think that 
this is due to r not being large enough. However, we found that deviations 
of f T (p) from p for large p are generic, while the straight line is approached 
only for smaller p. Nonetheless, we cannot say anything conclusive about the 
large r, large p behavior because the distribution gets very peaked and hence 
becomes smaller for large deviations, which makes the numerical method 
unreliable. 
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Figure 1: Numerically obtained f T (p) (bold dashed) for v* = 1.5 and r = 1.3, 
versus the (v* independent) extension of the SSFT for r — ► oo (bold solid). 
Also plotted are the conventional SSFT (thin solid), and the numerically 
obtained distribution function P T (pwr), scaled by its value at zero (thin 
dashed) . 

3.2 Large deviation method — extended SSFT 

Therefore, we used next an analytical asymptotic approach based on large 
deviation theory jH] similar to the treatment by Lebowitz and Spohn[5 a . One 
considers then 

e(A) = lim — -\n(e~ XQT ). (12) 

t^oo WT 

This infinite-r quantity is used to reconstruct the distribution function of Q T 
for large r by setting 

P t (Qt) ~ exp [-wre(Q T /wr)\ , (13) 
where e(p) is the Legendre transform of e(A): 

e(p) = max[e(A) — \p\. (14) 

A 
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conventional SSFT 

finite v. f x (p) ■ 

P z (pwx)/P x (0) " 



For a class of models, Lebowitz and Spohn proved the symmetry relation 



e(A) = e(l-A). (15) 

From this, using equations (JTTJ), (JT3|) and (JT3J), one sees that lim T _ )00 f T (p) = 
p, i.e., the conventional SSFT holdsjSJ. 

Our numerical results suggest, however, that for our model the conven- 
tional SSFT for the heat does not hold. We therefore expect equation ([To]) 
to be violated. Indeed, the following calculation of e(A) shows this to be the 
case. 

The Fourier transform P T (q) in equation (fTUj) determines e(A). First, from 
equation (JHJ), we have 

/oo 
dQe- XQ ^P T (Q T ) =P T (i\), (16) 
-oo 

Thus we need the analytic continuation of P T to imaginary arguments. This 
poses no difficulty as long as P T remains analytic. One finds from equa- 
tions (HH and (dnj) 



^- XQt ) = ~ 77 z " — ■ ( i7 ) 



exp 



«A(1-A)fr+^^ 



[1 - (1 - e~ 2 -)A 2 ] 3/2 

Clearly, there are divergences at the singular points A = ±(1 — e _2r ) -1 / 2 , 
where the right-hand side (r.h.s.) of equation (fTB*j) is no longer analytic, 
so that equation ([T7|) only holds for values of A in between those. Using 
equations (fT2|) and (fTTj) . we have 

e(A) = A(l - A) for |A| < 1, (18) 

where, taking r — ► oo as in equation ()12|) moves the singularities to ±1. This 
e(A) satisfies equation (fTB]) for < A < 1. 

However, as A approaches the singularities, the function in equation (fTTj) 
diverges. Beyond the singularities at ±(1 — e~ 2r )~ 1 / 2 , the r.h.s. of equa- 
tion (fT7j) becomes purely imaginary, and multivalued due to the denomi- 
nator. But the left-hand side of equation ()17|) remains real. Clearly, we 
cannot use equation (fT7|) for |A| > (1 — e _2r )~ 1//2 . To determine (e _A< ^) in 
that case, we first need to know why the integral in equation ()16|) diverges 
as A — > ±(1 — e _2r ) -1 / 2 . As we will argue next, this happens because P T 



7 



has exponential tails. Since P T (Q T ) is a normalized distribution and e~ x ® T 
a regular function, the divergence in equation (|TT)j) can only be due to the 
behavior of the integrand at ±00. In fact, for A > 0, any divergence must 
be due to the behavior at negative Q T and for A < 0, it must be due to the 
behavior at positive Q T . Now, for A < 0, if the distribution function P T (Q T ) 
fell off faster than exponential for large positive Q T , the factor e~ A< ^ T could 
not make the integral diverge. But as it does diverge, we conclude that the 
distribution function falls off exponentially or slower. On the other hand, 
if it did fall off slower than exponential, then the exponential factor e~ x ® T 
would always dominate the distribution function for large positive Q T and 
the integral in equation ()16)) would diverge for all A < 0. Since there are 
negative values of A for which the integral converges, the function P T cannot 
fall off slower than exponential. Hence, it must fall off exponential for large 
Q T . Considering A > 0, one deduces along similar lines that it must fall off 
exponentially for large negative values of Q T as well. 

In fact, the integral in (JT5j) diverges for all |A| > (1 — e~ 2r )~ 1 ^ 2 . For if the 
function P T (Q T ) falls off exponentially for large positive Q T , say as e _aQr , 
the integral in diverges for all A < —a. Likewise, given that the P T (Q T ) 
falls off like e a ® T for large negative Q T , the integral diverges for all A > a. 
Hence, for |A| > (1 — e _2r ) -1 / 2 , the quantity on the r.h.s of equation (fT2"|) 
of which the limit is taken, is minus infinity for all r, so that e(A) = —00. 
Thus, equation (|P%|) becomes 

e(A) = ( A ^- A ) f( ^J A ' < 1 (19) 



-00 otherwise. 

This e(A) does not satisfy the symmetry relation in equation 1)150 . e.g., for 
A = —1/2, e(A) = —3/4 whereas e(l — A) = —00. The fact that equation (JT5j) 
is not satisfied, means that the conventional SSFT does not hold. To make 
this more precise, we use equations (fT4*j) and (JTHJ), to find 

-p for p < — 1 

e(p) = { (p - l) 2 /4 for -1 < p < 3 (20) 
p — 2 for p > 3. 



•1 



Note that via equation (|T5j) . the large \p\ behavior is indeed exponential 

4 Singularities affected the distribution of a different quantity as well in J. Farago, J. 
Stat. Phys. 107, 781 (2002). 
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Using equations (jTTj) . (|T5|) and (J2UJ) . we find 



for < p < 1 

lim / T (p) = { p - (p - l) 2 /4 for 1 < p < 3 (21) 

for p > 3. 

For negative p, we have 

Zr(-P) = ~fr(p). (21') 

Equation (|21|) is an extension of the conventional SSFT. It coincides with 
it for the middle region —1 < p < l, 5 but differs from it for other p values. 
Most notably, for p > 3, it attains a constant value of two. 



3.3 Saddle-point method — large finite times 

If we compare the exact prediction of equation (J2TJ) (plotted as the bold 
solid line) with the numerical results (bold dashed line) in figure ^ a clear 
discrepancy emerges: The curve of f T keeps increasing with increasing p, 
whereas equation (}2*T]) predicts that it should level off to a value of two. This 
turns out to be a finite r effect. To prove this, we need a better treatment for 
large but not infinite r. This can be obtained from a saddle-point method 
applied to e(A), which we will present in a future publication The saddle- 
point method gives reliable results for sufficiently large r, as can be verified 
by a comparison to our numerical results^]. The asymptotic behavior for 
large r is then given by 

fW _/p+ %)/r + Q(r - 2 ) for p < 1 

1 ' I 2 + v/8(p-3)/r + Oir- 1 ) for p > 3, 1 ] 



where 



h(p) = — — In 

9 — p 2 2w 



(3-p)(l+p) 



.(3 + p)(l-p). 

We left out the behavior in between p = 1 and p = 3 as it does not appear to 
give any additional insight. Equation (}2*2*|) shows that as a function of p, the 
function f T (p) increases like ~ ^p — 3 (for fixed r), while as a function of r, 
it decreases (for fixed p). As expected, for r — > oo, it approaches the large 
deviation result in equation (J21j) as 1/r for small p and as 1/ y/r for large p. 

5 A restriction to the SSFT (i.e., \p\ < 1) was also found in a model of heat conduction 
in reference 
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Whether the new features beyond p = 1 are observable depends on the 
values of P T (pwr) and P T (—pwr). If these are too small, the corresponding 
f T (p) cannot be seen in an experiment. But the value of the distribution 
function as plotted in figure [T] is non-negligible for values for p (and —p) at 
which f T bends away from the conventional SSFT. Furthermore, the values 
v* = 1.5 and r = 1.3 used in figure Q are realistic, as in the experiment of 
Wang et al. v* ~ 2.5 and r goes up to ~ 6. So this behavior should be 
experimentally detectable. 

4 Discussion 

Summarizing, we have shown that the behavior of heat fluctuations in a dis- 
sipative system with a deterministic, mechanical component (the potential), 
and a stochastic component (the heat bath, i.e., the water), differs from that 
known from both purely deterministic systems and purely stochastic systems, 
in two respects, a) For infinite r, the behavior of the conventional SSFT is 
only seen for the scaled heat fluctuation p between —1 and 1. For p > 1, after 
a parabolic region between p = 1 and 3, the quantity f T no longer increases, 
but stays at a plateau value of two [for p < — 1, see equation (|21'j) ]. b) The 
finite t behavior of the conventional SSFT is in general unknown, but in our 
case we find that f T keeps increasing with p. However f T {p) stays well below 
the conventional SSFT , implying a larger ratio of the probabilities of the 
particle to absorb rather than supply heat. These features are observable. 

One of the striking features of the extended SSFT, the plateau value of two 
for large (infinite) r and large p, can be understood physically. For large r and 
large Q T , i.e. p, the exponentially distributed AU T far outweighs the Gaussian 
distributed W T in equation (0). The distribution of AU T is exponential for 
large values (ex e~^ Au ^), because it is the difference of the potential energies 
of the particle at two times [cf. equation (J2J)], which are both Boltzmann-like 
distributed (because of the presence of the water) and independent of each 
other for large r. As (Q T ) = wt, this leads to P T (Q T ) oc e~^ T ~ WT \ which 
yields P T (Q T ) / P T (-Q T ) w e 2f3wT (if Q T > 0). Then equation JTIJ) and (3 = 1 
give f T ~ 2 of equation (|21|). 

We only considered an extension of the SSFT here. An extension of the 
TFT can also be obtained. While for r — > oo, one gets again equation (|2*T|). 
for finite times, the extended TFT and SSFT differ |T3j. Furthermore, the 
extended TFT differs fundamentally from the conventional TFT, in that the 



10 



latter holds as an identity for all r, whereas the former one holds only in the 
t — > oo limit 6 . 

Finally, one may wonder about the generality of the extended SSFT. Our 
Langevin-based theory is only applicable near equilibrium. The arguments 
given above suggest that the extended SSFT could hold for other potentials 
as well. Perhaps it could even hold for a larger class of systems, not near 
equilibrium, with mechanical and stochastic components, since these are the 
main physical ingredients in our theory. 
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